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ABSTRACT 

We present the results of far-infrared imaging of extended regions around 
^ _ three bipolar pre-planetary nebulae, AFGL 2688, OH 231.8+4.2, and IRAS 

■ 16342—3814, at 70 and 160 fim with the MIPS instrument on the Spitzer Space 
^ ■ Telescope. After a careful subtraction of the point spread function of the central 

star from these images, we place constraints on the existence of extended shells 
and thus on the mass outflow rates as a function of radial distance from these 
QQ ■ stars. We find no apparent extended emission in AFGL 2688 and OH 231.8+4.2 

^ ' beyond 100 arcseconds from the central source. In the case of AFGL 2688, this 

O ' result is inconsistent with a previous report of two extended dust shells made on 

■ the basis of ISO observations. We derive an upper limit of 2.1 x 10^^ Mq yr~^ 
^ ■ and 1.0 X 10"^ Mq yr'^ for the dust mass loss rate of AFGL 2688 and OH 231.8, 
O ■ respectively, at 200 arcseconds from each source. In contrast to these two sources, 

IRAS 16342—3814 does show extended emission at both wavelengths, which can 
^ ■ be interpreted as a very large dust shell with a radius of ~ 400 arcseconds and 

■ a thickness of ~ 100 arcseconds, corresponding to 4 pc and 1 pc, respectively, 
at a distance of 2 kpc. However, this enhanced emission may also be galactic 
cirrus; better azimuthal coverage is necessary for confirmation of a shell. If the 
extended emission is a shell, it can be modeled, with some assumptions about 
its dust properties, as enhanced mass outflow at a dust mass outflow rate of 
1.5 X 10"^ M0 yr~^ superimposed on a steady outflow with a dust mass outflow 
rate of 1.5 x 10~^ M0 yr~^. Because of the size of the possible shell, it is likely 
that this shell has swept up a substantial mass of interstellar gas during its ex- 
pansion, so these estimates are upper limits to the stellar mass loss rate. We flnd 
a constant color temperature of 32 K throughout the circumstellar envelope of 
IRAS 16342—3814, which is consistent with heating by the interstellar radiation 
fleld. 
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Suhject headings: (ISM:) planetary nebulae: individual (AFGL 2688, OH 231.8+4.2, 
IRAS 16342-3814) — stars: AGB and post-AGB — stars: mass loss 



Introduction 



Pre-planetary nebulae (PPNe) represent a fleeting stage of stellar evolution. These 
transitional objects arise between the rapid mass loss phase at the end of the asymptotic 
giant branch (AGB) and the ionized planetary nebula stage. AGB stars are important 
because they are the Galaxy's main mechanism for the replenishment of dust and gas into 
the interstellar medium (ISM), ejecting most of their main sequence mass within a few times 
10^ years because of their high mass loss rates. By looking at the circumstellar material 
around PPNe, we can see an imprint of the activity of the earlier AGB. PPNe are also in 
a stage where the geometry of the mass loss from the central star usually changes from 
spherically symmetric to an axially sym metric bipolar outflow, as illustrated in the HST 
images of AFGL 2688 (ISahai et al.lll998l ). In these images, a spectacular bipolar structure is 
superimposed on a set of concentric shells of gas and dust ejected near the end of the AGB. 
In the study reported here, we have attempted to characterize the mass loss behavior of the 
progenitor AGB stages of three systems that are now bipolar. 

Using the direct detection of infrared emission from dust lost over as much as 10^ years 
during the AGB, it is possible to study the mass loss histories of these objects. There has 
been a dearth of observations of large extended shells of evolved stars out to the distances 
sampled in this project primarily because of the difficulty in observing in th e far infrared 



wher e the cool dust has its peak emission. With the Spitzer Space Telescope ( IWerner et al. 



20041 ). the far infrared is now accessible at a resolution and sensitivity sufficient to potentially 
resolve the structures of these dust shells, and thereby to test theories of mass loss in AGB 
stars. 

Several tantalizing cases were observe d by IRAS, showin g that large dust shells exist 
around some evolved stars. For example, iGillett et al. I Jl986h found that R CrB shows a 
very large shell with a radius of 4 parsecs, and iHawkind (Il990l ) reported a dust shell with 
a diameter of 30 - 40 arcminutes (~ 1 pc) around W Hya. However, these observations 
could not resolve the structure that may be present in the shells, or address the possibility 
of multiple shells. If the mass loss during the AGB is constant, we should expect to see a 
smooth envelope with column density declining as 1/6, where h is the displacement from the 
star in the sky. However, if the mass loss rate fluctua tes as the star goes throug h thermal 



pulses during the AGB, as the models suggest (e.g., IVassiliadis fc WoodI Il993l ). then we 



should see enhanced emission in the form of multiple dust shells around our objects. For 
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example, ISpeck et al.l (120001 ) reported on the basis of ISO data that there are large concentric 
shells in AFGL 2688 and AFGL 618. This result is among a small set of observations 
showing evidence for periodic mass loss caused by thermal pulsation as well as the time scale 
between thermal pulses predicted by AGB evolutionary models (e.g., IVassiliadis fc Wood 



19931 ). Confirmation of these shells at higher resolution and greater sensitivities motivated 
our selection of AFGL 2688 as a target for this study. There have been no previous claims 
of large-scale dust emission from OH 231.8+4.2 (hereafter OH 231.8) or IRAS 16342-3814. 
However, they are presently very young PPNe with high mass loss rates of approximately 
1O~^M0 yr~^ (lAlcolea et al.ll200ll : ISahai et al.lll999l ). which suggests that the emission from 
dust produced during the AGB may be readily visible. 

Here we report on the lack of very extended emission, spherically symmetric or otherwise, 
from AFGL 2688 and OH 231.8. However, we do see possible evidence for a very large diffuse 
dust shell around IRAS 16342—3814 with a radius of 400 arcseconds, though this extended 
emission might also be from galactic cirrus. 



2. Observations and Data Reduction 

AFGL 2688 ( the Egg Nebula), OH 231.8, and IRAS 16342-3814 were observed with the 



MIPS instrument ( iRieke et al.ll2004l ) at 70 and 160 /im on the Spitzer Space Telescope. These 
PPNe were observed along two mutually perpendicular scan paths in order to sample dust 
emission out to about 800 arcseconds from the central source, and to determine background 
levels. At 70 /im, the scan paths are 15' x 3.0' in one direction and 11' x 7.7' in the other, while 
at 160 fim, the scan paths are 15' x 2.6' and 10.5' x 6.5'. The pixel scale is 9.2 arcseconds 
at 70 /xm and 16 arcseconds at 160 /im. However, the central sources were not directly 
observed in order to avoid persistence artifacts from these bright sources. This observing 
strategy limits the azimuthal coverage of the sources, but should give an adequate estimate 
of the presense of well-defined shells or asymmetries such as large-scale bipolarity. 

The basic science calibrated data (BCD) was reduced using the software package MOPE}(0 
from the Spitzer Science Center (SSC). MOPEX was used for rejecting outliers as well as 
for co-adding and mosaicking the individual BCD frames to create the final image. Further 
reduction was done manually to remove striping due to bright latents as described in the 
MIPS data handboolj^. The correction involved finding the median value measured for every 



^http:/ /ssc. spitzer. caltech.edu/postbcd/ 
^http://ssc. spitzer. caltech.edu/mips/dh/ 
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pixel from a series of BCD frames far from the central source, where we assume there to be 
no extended emission, then subtracting this median from each data frame before co-adding 
and mosaicking. The median subtraction removes both the bright latents and the uniform 
component of background emission. Figure [T] shows the final mosaicked images. 

The most prominent features of the 70 fim images of the bipolar targets are the diffrac- 
tion spikes from the point spread function (PSF) of the central source. Although we did 
not image the central star, the wings of the PSF are still present and bright out to about 
200 arcseconds from the source. At 160 fim, the PSF is less pronounced but we see much 
more background than at 70 /xm from galactic cirrus. The galactic cirrus emission presents 
a problem in determining the true sky background since it covers much of our fields. For 
both the 70 and 160 fim images, we use the median surface brightness of regions of low 
and uniform brightness far from the central source as an estimate of the sky background. 
This procedure removes a uniform background from the images, retaining possible enhanced 
emission from the source and galactic cirrus. The uncertainty in the background (see below) 
at 160 is higher than at 70 fim because it is more difficult to find a large patch of uniform 
surface brightness to estimate the sky background. 



2.1. Sensitivities 

The sensitivity of our data is estimated by examining the standard deviation of a patch 
of uniform background far from the source both before and after the median background 
subtraction method described above. For AFGL 2688, before background subtraction, the 
mean of a patch in the eastern scan path at 70 /im is 14.1 MJy Sr~^ with a standard deviation 
of 1.11 MJy Sr~^. Some of the variance in the background is likely from the effects of bright 
latents from the detector pixels, which are removed by the median background subtraction 
described above. After median subtraction, the mean background is at 0.6 MJy Sr~^ with 
a standard deviation of 0.92 MJy Sr^^. The mean surface brightness at 160 fim before 
background subtraction is 31.45 MJy Sr^^ with a standard deviation of 1.67 MJy Sr^^. 
After background subtraction, the mean is at 0.33 MJy Sr~^ with a standard deviation of 
1.51 MJy Sr~^. The sensitivity for OH 231.8 is better than for AFGL 2688, probably because 
of the lower sky background levels. The sensitivity for IRAS 16342—3814 is comparable to 
that of AFGL 2688 at 70 fim, but worse at 160 /im due to a higher sky background and 
greater presence of galactic cirrus. The mean background and sensitivities for all three 
sources at 70 and 160 /im are given in Tables [1] and [2l respectively. 
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2.2. PSF Subtraction 

Since tlie PSF pattern from tlie central source is so prominent at 70 /xm, we attempted 
to fit a model PSF to our images in order to subtract the bright central source in search 
of extended dust emission near the source. The stinytim software provided by the SSC was 
used to generate model PSFs. We chose to use the default MIPS throughput curve along 
with a blackbody spectrum at several temperatures ranging from 20 to 100 K. See Figure 
|2]for examples of the model PSFs. Multidimensional fitting was carried out for our images 
in order to minimize the residuals of the PSF subtraction using various x-y shifts as well as 
the scaling of the PSF fiux values. The best fit was found by minimizing the square root of 
the sum of the squares of the differences in flux. We found several limitations in the model 
PSFs when attempting to subtract the PSF from the images. Even for the best-flt PSFs, the 
subtraction still leaves some obvious PSF residuals at a level of 10~^ relative to the peak as 
extrapolated from the model PSF. These residuals show that the wings of the model PSF are 
not well characterized at that flux level. The PSF subtraction also tends to leave portions 
of the region within the Airy ring with negative flux values. The most likely reason is a 
non-linear pixel response as the pixels approach saturation so that although the model has a 
higher brightness value closer to the core of the PSF, the actual pixel values are leveling off 
for our bright sources. Unfortunately, there are presently no empirical data on the MIPS 70 
or 160 /im PSF to compare to the model PSF at distances greater than 100 arseconds from 
the peak. Beyond about 240 arseconds, PSF features are no longer present in the MIPS 
images so our ability to detect extended emission there is only limited by background and 
integration time. 

The results of PSF subtraction for AFGL 2688 are shown in Figure|31 Note the remaining 
features near the north and at a position angle of ~ 50° east of north, about 150 arcseconds 
from location of the source in the PSF subtraction in Figure [3^. This residual can also been 
seen in the roughly the same region of the 70 /im PSF subtraction of OH 231.8 in Figure |H 
It seems unlikely that these features are physical, but the northern feature in AFGL 2688 
does align with its bipolar outflows seen in the optical and infrared, but which are observed 
within the inner 10 arcseconds of the PPN. The blob in AFGL 2688 at a position angle 
of ~ 50° also rough ly aligns with the outflow direction seen in previous CO observations 



( ISkinner et al.l 119971 ). also at around 10 arcseconds from the central star. The fact that the 
features we see are at nearly 200 arcseconds from the source and lie almost on top of the 
diffraction rays of the original PSF casts doubt upon their reality. This dilemma cannot 
be resolved until a well sampled empirical PSF out to several arcminutes is available. We 
may be seeing possible features in the wings of the actual PSF that the model has failed to 
reproduce. 
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As a substitute for an empirical PSF, we used OH 231.8 to subtract the PSF from AFGL 
2688. Figure [3)d shows the resuh. The residuals from using OH 231.8 as an empirical PSF 
appear to be much smaller than those obtained using the model PSF, reducing the residual 
of 30 MJy Sr~^ to less than 10 MJy Sr~^ at 100 arcseconds. Note that there are almost 
no PSF features, such as diffraction spikes, remaining in Figure [3)d compared to Figure [3^1. 
However, it may be problematic to use OH 231.8 +4.2 as an empirical PSF because it may 
have extended emission as well, although it would be remarkably fortuitous if both sources 
had extended emission having precisely the same morphology. It is also problematic to scale 
OH 231.8 to the same surface brightness as AFGL 2688 because in order to subtract the PSF, 
the whole image needs to be scaled, thus scaling the background as well. Because AFGL 
2688 is over twice as bright as OH 231.8, but the background in AFGL 2688 is not, scaling to 
the level of the AFGL 2688 PSF would overemphasize the background in OH 231.8 so that 
the residuals after the PSF subtraction might be dominated by the background in OH 231.8. 
We can nevertheless see that the PSF subtraction using OH 231.8 results in residuals at a 
level ~ 10 times below that from using the model PSF at 150 arcseconds, which suggests 
that OH 231.8 is either a point source at 70 /im or that its extended structure has the same 
orientation, scale, and shape as in AFGL 2688. 

The two-dimensional PSF subtraction at 160 /xm was done in a similar way using the 
model stinytim PSF. There are PSF features remaining after this subtraction as well (see 
Figure [3]); the ends of the diffraction spikes at 200 arcseconds from the source in the original 
image are not completely removed. We did not attempt a PSF subtraction of the 160 /im 
image of IRAS 16342—3814, because that the PSF of the central source is weak enough at 
this wavelength that its features are below the background and the extended structures. 

2.3. 1-D PSF Subtraction 

In order to investigate the hypothesis of spherical shells around a central source, we 
have determined an azimuthal average of the intensity around the central source. Such an 
average would show enhancements for spherically concentric shells projected as circularly 
symmetric features from episodic mass loss during the AGB phase. An azimuthal average 
also has the advantage that it greatly reduces the parameter space necessary to fit the model 
PSF. The source position and intensity are found by fitting to the brightest Airy ring visible 
in our images, which is the second brightest Airy ring of the PSF (see Figure [2]), since the 
source itself does not appear in our images. The PSF Airy ring used for scaling the PSF 
is located 76 arcseconds from the source at 70 fim and 170 arcseconds from the source at 
160 nm. We also find that the azimuthal averages provide a subtraction with less residuals 
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than using the two dimensional subtraction because localized irregularities of the PSF are 
averaged away. 

The one dimensional PSF subtractions for AFGL 2688 and IRAS 16342-3814 at 70 
fim were done using both the model PSF and OH 231.8, and are shown in Figure [51 The 
profiles of AFGL 2688 and OH 231.8 follow each other almost exactly out to beyond 200 
arcseconds from the source. This further indicates that OH 231.8 and AFGL 2688 are both 
strongly centrally concentrated at 70 /im, unless they both have exactly the same azimuthal 
excesses. They both also match the model PSF very well out to 100 arcseconds beyond 
which both appear to have slight excess (< 1% of the extrapolated peak of the PSF), above 
the model PSF. Once again, this excess in the wings is more likely due to the model PSF 
not reproducing the wings of the actual PSFs rather than resulting from physical emission 
associated with the sources. At 160 /xm, we only use the model PSF for subtraction because 
of possible contamination due to galactic cirrus in OH 231.8. 

3. Results and Analysis 
3.1. AFGL 2688 



3.1.1. 70 fim 

At 70 fim, the PSF-subtracted images of both AFGL 2688 and OH 231.8 are similar. The 
PSF residuals extend to about 150 arcseconds from the source, which may be obscuring some 
physical extended features. However, there are no obvious circularly symmetric residuals 
in the PSF subtraction even with the confusion near the central source, which would be 
expected if there had been spherically symmetric mass loss above a dust mass loss rate of 



2.1 X 10 ^Mq yr ^ (see discussion section). Previous studies with ISO by ISpeck et al. 



( I2OOOI ) led to the claim of large shells around AFGL 2688 having radii of 150 and 300 
arcseconds on the basis of a one-dimensional scan of the Egg Nebula with ISOPHOT at 
120 and 180 /im. We find, with better spatial resolution and azimuthal coverage using 
Spizer, that at 70 /xm, there are no signs of shell-like extended emission in our field. At 150 
arcseconds from the source, the PSF subtraction residuals have a surface brightness of about 
16 MJy Sr^^ using the model PSF and about 1 MJy Sr~^ using OH 231.8 as the empirical 
point source(Figure E]). Note also that there is no discernible excess emission above the 
fluctuations in the background (cr ~ 0.92 MJy Sr~^) at 300 arcseconds from the source. 
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3.1.2. 160 

The 160 /iin image (see Figure [3]) shows roughly uniform enhanced emission to the edge 
of our field of view at 1000 arcseconds from the source in the eastern scan path; the northern 
scan show enhanced emission out to about 300 arcseconds, then dropping off to background 
levels. The enhanced emission in the two paths is asymmetric as it does not fall off with 
distance in the eastern scan path. There does not appear to be any structure to the emission 
in either scan, other than that which may be attributable to the PSF diffraction spike in 
the north. Because the residual extended emission shows no symmetry about the central 
star, and because it shows no intensity falloff in the eastern scan path, this emission most 
likely arises from irregularly distributed galactic cirrus in this direction; this is supported by 
the relatively lower quality IRAS 100 /im image, which shows substantial extended cirrus 
beyound our survey region. The ISOPHOT observation of a possible shell at 300 arcseconds 
was based on a 53' x 3' linear scan with 30" x 92" pixels, at a position angle of 8 degrees east 
of north centered on the source, which in our case would be sampling the emission in the 
northern scan. We see with better azimuthal coverage with MIPS that the emission in both of 
our scans is likely strongly contaminated with galactic cirrus. The presence of the inner shell 
reported by Speck et al. at 150 arcseconds cannot be directly confirmed by our observations 
because the 160 /xm PSF Airy ring at 170 arcseconds from the source overlaps that region. 
The surface brightness at the Airy ring is about 17 MJy Sr~^ after background subtraction 
(but before the PSF subtraction). This is comparable to the backgrou nd subtracte d surface 



brightness of 20 (30) MJy Sr-^MJy Sr'^ at 120 (180) /xm reported bv ISpeck et al.l fl2000h . 



In the azimuthally averaged surface brightness of our 160 /xm image (Figure [6]), we see 
a similar result as the 70 /xm data. Speck et al. reported a brightness value of the putative 
shell at 300 arcseconds from the source to be 20 MJy Sr~^ at 120 yum and 60 MJy Sr~^ 
at 180 /xm, without background subtraction. The 120 /im emission is about 10 MJy Sr^^ 
above their background as extrapolated from the surface brightness far from the star in their 
plot. At 180 yum, the excess emission is about 15 MJy Sr^^ above the background. In the 
MIPS 160 /im image, we find that the surface brightness at 300 arcseconds from the source 
is about 2 MJy Sr~^ after background subtraction. Although the comparison is not at the 
same wavelength, the excess emission we detect is 7 times below that measured by ISO. 
When comparing the two scans in our observations separately, we see that the eastern scan 
has roughly uniform brightness of 4 MJy Sr~^ out to 1000 arcseconds, while the emission in 
the northern scan drops off at 400 arcseconds from the source from about 2 MJy Sr~^ to the 
level we use as the background in the data reduction. Both the asymmetry in the enhanced 
emission in both scans and the low surface brightness in the Spitzer data at 300 arcseconds 
from the source lead us to conclude that the shell at this distance reported on the basis of 
ISO data is probably not associated with the source. 
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3.2. OH 231.8+4.2 

Neither the 70 nor 160 fim images of OH 231.8 show any signs of spherically symmetric 
extended emission. The PSF subtraction at 70 fim leaves similar residuals as for AFGL 
2688, though the brightnesses of these residuals are much lower because the central source is 
intrinsically less bright. The strongest PSF residuals are along the diffraction spikes, similar 
in shape to those from AFGL 2688. This further suggests that these residuals are not physical 
features, but are rather PSF artifacts not accounted for by the stinytim model. In the 160 
fim azimuthally averaged plots, both scans appear to have rather uniform emission, with the 
northern scan having slightly higher surface brightness at approximately 2 MJy Sr~^, while 
the eastern scan has an average surface brightness of 1 MJy Sr~^ (Figure [7]). 

The 160 fim image shows some slightly clumpy emission along the northern scan path, 
slightly beyond the PSF diffraction spike (see Figured]) at 320 arcseconds. The clumps, with 
a width of about 80 arcseconds, have a surface brightness of about 3.6 MJy Sr~^ on a 2 
MJy Sr~^ diffuse background. The IRAS 100 /xm image shows diffuse galactic cirrus in this 
region, which may be the origin of these clumps. They are aligned with the bipolar axis, but 
unfortunately also with the diffraction ray so it may also be a PSF feature. Since this study 
is concerned with investigating spherically symmetric emission, we will defer the question 
of their nature. The azimuthally averaged surface brightness shows no enhanced emission 
attributable to the source (Figure [H]). 



3.3. IRAS 16342-3814 

In contrast to the previous two sources, IRAS 16342—3814 does show some evidence for 
what may be a large patchy shell at both 70 and 160 fim. The enhanced emission out to 
a radial distance of ~ 400 arcseconds is consistent with rough circular symmetry about the 
star in the coverage area available (Figured]). This radius corresponds t o a physical size o f 



about 4 parsecs, assuming the distance to IRAS 16342—3814 to be 2 kpc (ISahai et al.lll999l ). 
This possible shell structure is rather patchy, with a large arc in the western scan at 70 fim. 
At 160 /xm, the same arc is present, but is more diffuse. Nevertheless, the azimuthal average 
shows an excess above the background at ~ 400 arcseconds from the central source. The 
brightest portion of this feature occurs at a position angle of 250 degrees east of north, which 
is the same direction as the axis of the bipolar nebula seen in the optical. The ratio of the 
70 to 160 /xm average surface brightness, is roughly constant (~ 0.3) as a function of radial 
distance from the source. This ratio yields a color temperature of ~ 32 ± 2 K, and, if a z/^ 
emissivity law is assumed, a dust temperature of 26 ± 2 K. The constant color temperature 
suggests that the dust is heated primarily by the interstellar radiation field. 
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Because there are a number of interstellar cirrus features in the general direction of 
IRAS16342-3814, we cannot definitively argue that the patchy arc of emission in the scan 
paths around this source is caused by mass loss from the star. The IRAS 100 /im image of 
this region with an overlay of the coverage of our data (figure ^ shows that there is indeed 
a substantial likelihood that this emission is cirrus. With this caveat in mind, we proceed in 
the analysis below on the assumption that the extended emission has resulted from mass loss 
from the star to obtain some indicative numbers for the mass loss rate necessary to produce 
such a shell. In the following analysis, we will use the term "mass outflow rate" to indicate 
the total mass flowing outward per unit time across a spherical shell at a particular radius; 
this includes both the stellar component and the interstellar medium being swept up; we 
will use the term mass loss rate to refer to mass loss from only the star. 

Making a few assumptions about the mass outflow and the dust, we can set some limits 
on the mass outflow rate of IRAS 16342—3814. Assuming that the mass loss is spherically 
symmetric, we can fit a model of the expected brightness profile to the observed radial profile 
of the source. Figures [TO] and [TT] show the azimuthally averaged radial profiles at 70 and 160 
/^m, respectively. We fit only the 70 /xm radial profile because the re is less cont a minat ion by 



galactic cirrus compared to 160 fim. Following the derivation by iGillett et al.l (Il986l ). if we 
assume a 1/r^ density profile, appropriate for constant mass outflow, where r is the distance 
from the star, and a constant expansion velocity, Vg, then the surface brightness Iu{b) has 
the form: 



iu\b) = z ^ cos ' 



27CVpb 



R„ 



2TTVeb 



for R. 
b 



cos \- ) - cos \ 

^max Rrnin 



T, ^ b ^ Rjnax 
for b < Rrr,.ir. 



where b is the projected physical distance from the source in the sky, is the opacity 
coefficient in units of cm^ g~^, By{T) is the blackbody function, and Rmax and Rmin are the 
outer and inner radii of the shell. We further assume that the outflow velocity is constant 
at 15 km s~^, and that the dust is composed of small astronomical silicate grains with Ki, 
= 104 cm^ g~ ^ at 70 /xm (using a = 0.1 fim, p = 2.3 g cm~^, and QJ^^^ = 2.99 x 10~^ from 
Draine fc Led (119841 )). We also adopt a constant temperature of 26 K throughout the shell 
and envelope. Figure [T^] gives a comparison between different models having various mass 
outflow rates. The best-fit model is a shell with a 4.2 pc radius and a thickness of 1 pc, with 
a dust outflow rate of 1.5 x 10^^ Mq yr~^ superimposed on a smooth envelope with a constant 
dust outflow rate of 1.5 x yr^^. This implies a rather high gas outflow rate on the 

order of 3 x 10"'^ Mq yr~^ in the shell and 3 x 10~^Mq yr~^ in the smooth envelope, assuming 



- 11 - 



a gas to dust mass ratio of 200. Using OH 231.8 as the empirical PSF, the fit only requires a 
shell component with a mass outfiow rate of 3 x 1O~^M0 yr~^ without an underlying smooth 
envelope. These fits depend on the choice of PSF, but once the PSF is chosen, the mass 
outfiow rates are insensitive to residuals from PSF subtraction, since we only fit for the 
region beyond 150 arcseconds from the central source (see Figure [TOl) . However, our value 
for the mass outfiow rate is highly dependent on the assumed temperature of the dust. For 
typical ISM dust temperatures between 22 and 35 K, the inferred dust outfiow rate is related 
empirically to the adopted temperature by: M oc T^''-^ at 70 /xm. If T > 26 K, we would 
infer a substantially lower mass outfiow rate. 

Given the model thickness of the shell and assuming a typical expansion velocity of 15 
km s~^, we estimate the maximum duration of the enhanced mass loss event which produced 
the shell to be about 65,000 years. This is an upper limit because the internal velocity 
dispersion in the shell broadens the shell as it expands. A 1 km s~^ dispersion, for example, 
would reduce the duration of the mass loss that produced the shell by about a third. If we 
use an age of 40,000 years for the duration of the enhanced mass loss and approximate the 
shell as spherical, the total dust mass would be about 0.04 Mq, giving a total shell mass of ~ 
8 Mq, assuming a gas-to-dust mass ratio equal to 200. The velocity of the shell is probably 
slower than the molecular outfiow velocity of a typical AGB star and the total amount of 
mass lost by the AGB star is likely smaller than this amount because of the interaction with 
the ISM out to the distances we are observing (see below). 

The mass of the shell can be estimated directly from the infrared emission using the 
equation M = F^D'^ / {By{T)Ky), where is the total fiux density in the shell and D is the 
distance to the source. The integrated fiux density from the limited azimuthal coverage that 
we have between radii of 300 and 400 arcseconds is 6.4 Jy at 70 /im, corresponding to a dust 
mass of 0.015 in our observed portion of the shell, assuming a temperature of 26 K for 
the dust. We have about one fifth of the full azimuthal coverage with data at this distance 
from the star, and can estimate the total dust mass of the shell if we approximate the shell 
as isotropic with an averaged surface brightness of 2.5 MJy Sr~^ at 70 /im from the radial 
profile, with a width of 100 arcseconds, and a temperature of 26 K. The total dust mass of 
the extrapolated shell is then about 0.03 Mq, comparable to the above estimate. 

Because the shell is so large, the amount of interstellar matter that may have been swept 
up and accumulated in the shell during the AGB could be significant. For an ISM hydrogen 
density of 1 cm~^, a 4 parsec radius shell would have swept up about 7 Mq, which would 
account for almost all the mass that we may be measuring. The accumulation of interstellar 
material could potentially lower the stellar mass loss rate estimated from our model by a 
large factor. We therefore emphasize that the mass outfiow rate estimate is only an upper 
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limit to the stellar mass loss rate. 

One of the difficulties in addressing the mass loss history of IRAS 16342—3814 with the 
current data is its location in the Galaxy; with a scale height of only 150 pc, there is possible 
confusion with diffuse galactic cirrus. An unfortunate alignment with background emission 
is an alternative to the existence of a shell produced by mass loss. Further observations with 
better azimuthal coverage are needed to test the circumstellar shell hypothesis. 



4. Discussion 



The extended dust emission in AFGL 2688 and OH 231.8 seems to show that these two 
objects do not have as long a mass loss history as one might have anticipated. Our results 
for these two sources probe the region beyond ~ 100 arcseconds from the central star, which 
corresponds to 1.0 x 10^ yrs ago for A FGL 2688 an d 2.8 x 10'^ yrs ago fo r OH 231.8, given 
a distance of 420 pc (lUeta et al.ll2006l ) and 1.2 kpc (I Jura fc Morris! Il985l ). respectively, and 
an expansion velocity of 20 km s~^. We can can set a limit to the dust mass loss rate by 
inverting equation [TJ 

K^B^{T) COS ^{h/Rmax) 

The dust mass loss rate can be simplified as a function of the observed surface brightness I 
and the projected distance from the source in the sky h with the assumption that h <^ Rmax- 



70 urn 



M(70/im) ~ 3.8 x IO'^Mq yr'^ (^-^ 
M(160/im) ~ 2.9 x 10"^Mo yr'^ (^-^ 




MJy Sr 



30 K 



If 6 ~ Rmax, then the full form of equation [2] must be used because the cos~^ {b / Rmax) 
term in the denominator of equation [2] becomes important and will cause the dust mass 
loss rate inferred from a given surface brightness value to increase drastically. For the cases 
such as AFGL 2688 and OH 231.8 where we do not see a clear envelope associated with the 
source, we cannot be sure that the condition b -C Rmax holds since Rmax is indeterminate. 
The equations above also show the temperature dependence in the form of a power law to 
approximate the blackbody function between 22 and 35 K. From these equations we see that 
the temperature dependence at 160 /xm is less steep than at 70 /im, but the contamination 
by galactic cirrus near the sources makes estimating upper limits problematic at 160 fim. 
We can establish an upper limit to the mass loss rates for AFGL 2688 and OH 231.8 at 70 
fim based up on the surface brightness of the residual left from the PSF subtraction at 200 
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arcseconds from the source, which is a compromise between a location far enough from the 
source that PSF subtraction errors are small and being close enough to the source that there 
could plausibly be a circumstellar envelope. For AFGL 2688, 200 arcseconds corresponds to 
a radial distance f rom the source of 0.4 parsecs, using a distance to the source of 420 parsecs 
( lUeta et al.l l2006l ). The surface brightness after PSF subtraction with the model PSF is 
about 14 MJy Sr"^ a t 70 um. Assuming an expansion velocity of 20 km s ^ as observed in 
CO bv lSkinner et al.l (11997 ) , a du st temperature of 30 K, n^, = 104 cm^ g~^ appropriate for 
carbon dust (jPraine fc Led 1 19841 ). and that the possible envelope has a radius significantly 
greater than 0.4 pc, then the dust mass loss rate upper limit is 2.1 x 10^^ M0 yr^^. Similarly 
for OH 231.8, the residual emission from the PSF subtraction is 3 MJy Sr^^ at 200 arcseconds 
from the central star, which corresponds to 1.3 parsecs at a distance of 1.3 kpc. This surface 
brightness implies a dust mass loss rate of 1.0 x 10"^ M© yr ^ assuming an outflow velocity 
of 15 km s~\ Kj, = 98 cm^ g~^ (appropriate for silicate dust), and a temperature of 30 K. 

There are several potential explanations for the lack of extended emission seen in this 
study. One is that, even with Spitzer's increased sensitivity, the dust emission is below the 
threshold for detection. Since the distances we are studying are far from the central star, 
we would expect the temperature of the dust to be determined by the ambie nt interstellar 



radiation field. If the dust temperature is about 20 K, typical for the ISM (IMathis et al. 



19831 ). and an emissivity oc z/^, the surface brightness would be 12 times stronger at 160 
/xm, than at 70 fim where the resolution and sensitivity of MIPS is better. We would be 
more likely to detect cool dust emission at 160 fim than at 70 fim. However, confusion with 
the galactic cirrus is also likely at 160 fim because it is about the same temperature as any 
hypothetical extended dust emission associated with the source that is heated primarily by 
the interstellar radiation field. For comparison, our detection of the possible shell from IRAS 
16342-3184 has a ratio of Iwofj-m/hofim ~ 3.3, corresponding to a color temperature of 32 K, 
which requires a higher than average interstellar radiation field. 

The upper limits to the mass loss rates derived in this study can be compared with 
theoretical AGB evolutionary models, particularly those that predict enhanced mass loss 
rates due to thermal pulsation near the end of the AGB. The spatial coverage in this study 
between ~ 200 to 1000 arcseconds from the central star translates into a probe of the history 
of mass loss between 2 x 10^ yrs to 1 x 10^ years ago for AFGL 2688 and 6 x 10^ yrs to 
3 X 10^ yrs ago for OH 231.8, u s ing an expansion velocity of 20 km s~^ for both sources. The 
models by I Vassiliadis fc WoodI (119931 ) show that for a 2.0 M0 progenitor, during the last few 
xlO^ yrs of AGB evolution, there are several thermal pulses which result in enhanced mass 
loss rates peaking at about 1.3 x 10~^ Mq yr~^ during the end of each pulse. Using a gas to 
dust mass ratio of 200, we find that for AFGL 2688 and OH 231.8, the upper limit to the 
total mass loss rate during the above time intervals is about 4 x 10~^ Mq yr~^ and 2 x 10~^ 
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Mq yr~^, respectively. These limits are close to the sensitivity necessary to see shells that 
may be the result of the r mal p ulsation on the AGB. Although we do not detect the shells 
reported by ISpeck et al.l (I2OOOI ). which they attribute to thermal pulses, the signatures of 
thermal pulses may still be present, but below our current sensitivity. 

We also consider the possibility that OH 231.8 does not show extended emission in 
the far infrared because of interactions with a binary companion, which caused the central 
star to lose mass more rapidly and more recently than during the evolution of a lone AGB 
star. A companion to QX Pup (the central star of OH 231.8), is evidenced by optica l 
spectra consistent with a companion of stellar type AO V (ISanchez Contreras et al.l 120041 ). 
QX Pup also shows a peculiar paradox of being an M9III (jCohenlll98ll ) AGB star, while its 
bipolar activity and morphology display all the signs of post-AGB activity of typical PPNe. 
Having a close companion woul d enhance th e mass loss rate and provide a mechanism for 
generating the bipolar outflows (jMorrid 119871 ). thus shortening its mass loss history enough 
that we should not be surprised to se e no emission far fr om the source. In contrast to 
the colhmated outflow from OH 23 1.8 JAlcolea et al.l l200ll) . AFGL 2688 has a spherically 
symmetric envelope seen in ^^CO (lYamamura et al.l Il996) and evidenced by the partial, 
concentric, circular arcs present in HST images (ISahai et al.lll998l ). Since a binary interaction 
would not cause the past spherically symmetric mass loss, a possible companion is probably 
not a good explanation for initiating mass loss, although the present bipolar morphology of 
AFGL 2688 is consistent the possibility that a binary interaction has altered the mass loss 
in more recent times. 

Tracing the mass loss history during the AGB phase should be relatively straightforward 
via mapping the emission from the circumtellar envelope. But in practice, it has been 
difficult because molecular-line observations are ultimately limited by the photodissociation 
of molecules in the outer envelope regions, and far-infrared observations of dust emission have 
been limited by the generally low angular resolution of the space-based telescopes which have 
been available for this purpose in the past (IRAS and ISO). Hence the reported detections of 
very extended emission in a few dying stars with IRAS and ISO has generally been recognized 
as an important milestone in the study of mass loss. However, such detections have also raised 



thevery importa nt question of why only a few select objects like Y CVn (llzumiura et al. 



19961 ) or RY Dra (lYoung et al.lll993l ) - not particularly known for having high mass loss rates 
- show extended envelopes, whereas large numbers of stars having high CO-determined mass 
loss rates do not reveal the presence of such envelopes. Is this because the radial density 
law for most of these high mass-loss stars is significantly steeper than in objects like Y CVn 
or RY Dra, and in particular steeper than r^^ (an issue of profound importance for the 
evolutionary times of stars through the AGB phase and theories of mass-loss), or are the 
claimed detections of extended envelopes really a result of poorly characterised instrumental 
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artifacts? The Spitzer data presented in this paper clearly show that the presence of these 
shells cannot be confirmed at a level well below the intensities expected from Speck et al.'s 
results. Our non-detections call into question not only the ISO results on AFGL 2688, but 
all results on the detection of extended envelopes in other objects using the same linear scan 
technique described by Speck et al.. More detailed mapping of many high mass- loss rate 
objects like AFGL 2688, to search for such shells, is crucially needed. 



5. Summary 

Spitzer observations of extended envelopes of expanding, dusty outflows from bipolar 
AGB and post-AGB stars reveal that there may be a very large dust shell having a radius of 
400 arcseconds around IRAS 16342—3814. The combination of the presence of nearby cirrus 
emission and the limited azimuthal coverage in the images makes the conclusion for a shell 
uncertain, but if the shell is indeed the result of an episodic mass loss event, then it would 
represent one of the largest dust shells found so far. 

Our observations of AFGL 2688 at 70 um do not show the dust shells at 150 and 



300 arcseconds from the source reported by ISpeck et al.l (120001 ). Since the dust shell may 



be very cool, the non-detection at 70 /im does not alone rule out dust shells. However, 
we find that there is only a slight excess at 160 fim of 2 MJy Sr^^ above the background 
at the reported location of the outer shell (~ 300 arcseconds). With greater azimuthal 
coverage than was previously obtained with ISOPHOT we find that there is substantial 
contamination by galactic cirrus in the region at 160 /im, with galactic cirrus emission above 
the sky background present throughout the entire eastern scan path of our observations. 
Unfortunately, at 160 /zm we can only set an upper limit to the emission from a dust shell 
at 150 arcseconds from the source because one of the Airy rings of the PSF overlaps that 
region. We also see no extended emission from OH 231.8 at either 70 or 160 /im other than 
that attributable to galactic cirrus. In fact, using OH 231.8 as an empirical point source for 
PSF subtraction from AFGL 2688 at 70 fim appears to substantially reduce PSF subtraction 
residuals compared to using the model PSF. 

The limitation of our method of observation is that we have only two radial directions to 
probe possible extended emission. For cases like IRAS 16342—3814 where the shell is patchy, 
better azimuthal information would help to resolve whether the origin of the emission is from 
the star or from galactic cirrus. The MIPS data show that observing very extended emission 
in the far infrared is possible, but difficult because of the prominence of galactic cirrus 
emission at these wavelengths. 
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Mean Before 




Mean After 






Subtraction 


a 


Subtraction 


a 


Object 


(MJy Sr-i) 


(MJy Sr-i) 


(MJy Sr-i) 


(MJy Sr-i) 


AFGL 2688 


14.1 


1.1 


0.6 


0.9 


OH 231.8 


10.2 


0.9 


0.1 


0.6 


IRAS 16342-3814 


30.3 


1.3 


0.3 


0.9 



Table 1: The 70 /im mean surface brightness and sensitivity for a region of uniform back- 
ground far from the source, before and after the median background removal described in 
the text. 
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Mean Before 




Mean After 






Subtraction 


a 


Subtraction 


a 


Object 


(MJy Sr-i) 


(MJy Sr-i) 


(MJy Sr-i) 


(MJy Sr-i) 


AFGL 2688 


31.4 


1.7 


0.3 


1.5 


OH 231.8 


23.1 


0.7 


-0.2 


0.7 


IRAS 16342-3814 


90.1 


2.9 


1.0 


2.3 



Table 2: The 160 jim. mean surface brightness and sensitivity for a region of uniform back- 
ground far from the source, before and after the median background removal described in 
the text. 
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Fig. 1. — The three PPNe in our sample after mosaicking and background subtraction. The 
units of the colorbars are in MJy Sr~^. The bars in the images are 200" in length for scale 
and the x marks the location of the point source. Top: AFGL 2688. Center: OH 231.8. 
Bottom: IRAS 16342—3814. Note the prominence of the PSF features at 70 fim in all the 
images. At 160 /im, the PSF is not apparent in IRAS 16342—3814. North is up and east is 
to the left. 



Fig. 2. — Model stinytim PSFs in a log stretch to emphasize the Airy rings and diffraction 
spikes. Left: The 70 /im PSF. Right: The 160 /im PSF. The lines in the figures are 200 
arcseconds in length. For scale, the surface brightness of the diffraction spikes ~ 200 arc- 
seconds from the center is ~ 10~^ times lower than at the center. See Figures [3] and [5] for 
azimuthally averaged plots of the model PSF. 
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Fig. 3. — PSF subtracted images of AFGL 2688. The color stretches are hnear, with the 
units in MJy Sr^^. The x marks the location of the point source while the line represents 
200" for scale, a) PSF subtraction using the stinytim model PSF at 70 yum. The three PSF 
diffraction spikes visible in Figure [1] appear to subtract differently, leaving different levels 
of residuals, b) PSF subtraction at 70/im using OH 231.8 as a PSF, showing almost no 
residuals except for the region at about 60 degrees east of north, c) PSF subtraction using 
the stinytim PSF at 160 /xm. The feature in the northern scan path coincides with the 
position of the diffraction spike in the model PSF. The images have the same orientation as 
Figure [T], with north being up. 
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Fig. 4. — Left: PSF subtraction of OH 231.8 using the model PSF at 70 /im. The residuals 
are very similar to those resulting from the PSF subtraction of AFGL 2688 in Figure [31 
Right: PSF subtraction at 160 fim. The northern scan path has more structure than the 
eastern scan, probably from galactic cirrus. The units of the colorbars are in MJy Sr~^. The 
X marks the location of QX Pup, the central star of OH 231.8. North is up as in Figure [H 



-24- 
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Fig. 5. — Left: The azimuthally averaged surface brightness of AFGL 2688, OH 231.8, and 
the model PSF at 70 /im. The model PSF and OH 23L8 are scaled to the Airy ring from the 
AFGL 2688 PSF at 76 arcseconds. Right: The result of the PSF subtractions. The surface 
brightness drops to background levels beyond about 250 arcseconds. The excess emission 
left from the PSF subtraction using the stinytim model PSF is very similar to the residuals 
left after the same model PSF subtraction of OH 231.8 in Figure [71 The PSF subtraction 
using OH 231.8 as a template PSF instead of the model PSF is shown with a dashed line; the 
excess emission from the subtraction using the model PSF is almost completely eliminated. 
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Fig. 6. — Left: Azimuthally averaged surface brightness of AFGL 2688 at 160 /im. The 
two scan paths of the observations are plotted separately to show the differences in their 
background levels. The solid line represents the average of the two scans. Beyond about 
750 arcseconds from the source, the northern scan path stops so only the eastern path is 
represented beyond this distance. The PSF is scaled to the brightness of the Airy ring at 
170 arcseconds from the source. Right: The result of PSF subtraction. The region less 
than 200 arcseconds from the source may be affected by the PSF subtraction so we cannot 
make strong quantitative statements about the emission there at this time, except that the 
surface brig h tness of this region is likely to be about 2 MJy Sr^^. ISOPHOT data from 



Speck et al.l (120001 ) are also plotted for comparison after subtracting the background offset 



determined from the median of the values beyond about 400 arcseconds from the source. 
The reported ISOPHOT intensity at 300 arcseconds is almost eight times higher than the 
surface brightness seen by MIPS at 160 fim after PSF subtraction. 
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Fig. 7. — Left: the azimuthally averaged surface brightness of OH 231.8 and the model PSF 
at 70 /im. The model PSF has been scaled to match the Airy ring emission of OH 231.8 at a 
radius of 76 arcseconds. Right: the result of subtraction of the model PSF. The shape of the 
residuals is almost the same as those left from the PSF subtraction of AFGL 2688 (Figure 
|5]), falling off to background levels beyond 250 arcseconds from the source. 



OH 231 .8 1 60 microns Azimuthal Plots 

1 0000.0 M 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

; Source 

Eastern Scan Path 



0) 

c 



1000.0 r 



<n 100.0 r 



10.0 r 



1.0 



0.1 



Northern Scan Path 

Model PSF 




100 200 300 400 500 600 
Distance from Center (arcseconds) 



OH 231.8 160 microns PSF Subtracted 



10 



tn 

03 

c 



Average Azimuthal Brightness 

Eastern Scan Path 
Northern Scan Path 




100 200 300 400 500 600 
Distance from Center (Arcseconds) 



Fig. 8. — Left: the azimuthally averaged surface brightness of OH 231.8 along with the 
model PSF at 160 /im. Right: the 1-D PSF subtraction. There is no apparent enhanced 
emission other than attributable to fluctuations from the galactic cirrus. 
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Fig. 9.— The MIPS 70 //m field of view of IRAS 16342-3814 overlaid on the IRAS 100 
image of the surrounding regions. Note the abundance of galactic cirrus. The bar above the 
frame is 400 arcseconds in length, representing the radius of the potential shell seen with 
MIPS around IRAS 16342-3814. 
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Fig. 10.— Azimuthal average and PSF subtraction of IRAS 16342-3814 at 70 /im. Left: 
As with Figure O this plot shows the source plotted along with the radial profiles of OH 
231.8 and the model PSF. Unlike AFGL 2688, IRAS 16342-3814 shows a definite excess 
in its radial profile above the background and PSF profile. Right: the difference between 
using the model and OH 231.8 for the PSF subtraction. Both give about the same level of 
excess emission of ~ 2.5 MJy Sr^^ between 300 to 400 arcseconds from the source, where we 
see a patchy shell in the MIPS 70 and 160 /xm images. 
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Fig. 11. — Left: The azimuthally averaged surface brightness of IRAS 16342—3814 and the 
model PSF scaled to the brightness of the source at the location of the 170 arcsecond radius 
Airy ring of the 160 fim PSF. Right: the result of PSF subtraction. The excess emission 
seen in the MIPS images is clearly present out to 400 arcseconds. 
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Fig. 12. — Plot of the expected surface brightness for various mass loss rates. The solid line 
shows the expected surface brightness for a smooth constant mass loss envelope with a radius 
of 4.2 pc while the dashed lines show the radial profile of a single shell of 1 pc thickness at 
the same mass loss rate. The diamonds show the radial profile of IRAS 16342—3814 after 
the subtraction of the model PSF. We find that the best fit to the radial brightness profile 
obtained from subtracting the model PSF is a two-component model with a shell of radius 
4.2 pc and thickness 1 pc requiring a dust mass loss rate of 10~^ Mq yr~^ superimposed 
upon an envelope from a constant dust mass loss of 10~^ Mq yr^^ (dashed blue line). The 
asterisks show the radial profile found from using OH 231.8 as a substitute for the model 
PSF in the PSF subtraction will fit using only the shell component with a dust mass loss 
rate of 10~^ Mq yr~^ (dashed orange line). Points closer than a radial distance of ~ 150 
arcseconds are unreliable because they are most affected by PSF subtraction residuals and 
are not used in the fit. 



